Wetting regimes and interactions of parallel plane surfaces in a polar liquid 
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We apply a phenomenological theory of polar liquids to calculate the interaction energy between 
two plane surfaces at nm-distances. We show that depending on the properties of the surface-liquid 
interfaces, the interacting surfaces induce polarization of the liquid in different ways. We find, in 
full agreement with available experiments, that if the interfaces are mostly hydrophobic, then the 
interaction is attractive and relatively long-ranged (interaction decay length A ~ 1.2 nm). The water 
molecules are net polarized parallel to the surfaces in this case. If the surfaces are mostly hydrophilic, 
then the molecules are polarized against the surfaces, and the interaction becomes repulsive, but 
at a short-range (A ~ 0.2 nm). Finally, we predict there exists an intermediate regime, where the 
surfaces fail to order the water molecules, the interaction becomes much weaker, attractive and, at 
relatively small distances, decays with the inverse square of the distance between the surfaces. 



Interaction forces between hydrated, nm-size objects 
at short distances play an important role in various bio- 
logical and nano-fabrication processes. For example, the 
disjointing pressure between two biological membranes in 
pure water at distances 0.5 nm < d < 2.5 nm corresponds 
to a short-range repulsive force 



P 



-S- l dG/dd = P exp(-d/A), 



(1) 



where S is the cross-sectional area, G(d) is the inter- 
action energy of the system, 5 • 10 9 dyn/cm 2 < Pq ^5 
5 ■ 10 10 dyn/cm 2 , and A ~ 0.2 nm [I]. On the other 
hand, two hydrophobic plane surfaces exhibit attraction, 
P = — Pq exp (— d/X'), characterized by a similar prefac- 
tor, Pq ~ Pq, but different exponent: A' ~ 1.2 nm [2j[3]. 
Therefore, experiments show that hydration forces are 
characterized by at least two different length scales and 
depend on surface material properties. 

The nature of these forces are explained by a few theo- 
retical approaches. The Landau-type model with the or- 
der parameter corresponding to the ordering of the water 
molecules was presented in [3] to describe the repulsion of 
hydrophilic surfaces. In a related approach, this order pa- 
rameter was instead associated with the hydrogen-bond 
network deformations [5]. In [MS], the hydration forces 
between phospholipid membranes were associated with 
non-local polarization of the liquid. All of the models 
are purely phenomenological and provide an explanation 
for the repulsion force (fTl), although the specific values 
for the parameters Pi and A cannot be established from 
the theory. 

In this letter, we use a previously developed phe- 
nomenological model [9 -- 11 of a polar liquid to describe 
the interaction between plane surfaces arbitrary interface 
properties. We characterize the liquid by the microscopic 
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average of the molecular dipole moment orientations vec- 
tor, s(r), over a microscopic volume element centered at 
position r that contains a macroscopically large number 
of molecules. In the following, we represent the free en- 
ergy of the liquid as G[s (r)] = Gb + Gs, where Gs is the 
energy of the liquid-surface interface (see below) , and 



G B [s(r)]=P 2 JdV 



a,/3— x,y : z 



ds a ds a 

dx/3 dxp 



dV 



P 2 V(s 2 



8tt 



-El 



E e (r) 



(2) 



is the energy of the liquid bulk. Here, Pq = nodo, no is 
the density of the liquid, do is the molecular dipole mo- 
mentum and the Oseen energy term is characterized by 
C fw 0.5 nm 2 , the phenomenological constant responsible 
for the short range hydrogen bonds stiffness. The polar- 
ization density of the liquid, P = P s is related to the 
density of the polarization charges, pp = — divP. The 
polarization electric field is the solution of the Poisson 
equation divEp = — 47rpp, and E e is the external elec- 
tric field in the absence of the liquid. 

When the liquid polarization is small, s«l, "the equa- 
tion of state" function takes the usual Ginzburg-Landau 
form T^(s 2 ) « As 2 /2 + Bs A , where A and B are the phe- 
nomenological liquid-dependent constants. The former 
is determined by the long-range dipole-dipole interac- 
tions in the liquid and is related to dielectric constant, 
A = An/ (e — 1). Water is characterized by a large value 
of e sa 80 ^> 1. Therefore, A w 0.16, depends on the 
temperature and includes the entropy contribution aris- 
ing due to the averaging over the molecular orientation. 
The smallness of A is related to proximity of ferro-electric 
phase, A m (T — Tq) /Tq, predicted by the model and 
recently observed at temperatures T — Tc ~ 228K 
|101 [P2] , On the contrary, the parameter B depends 
on the short-range physics only, and is practically tem- 
perature independent, B ~ 1. The thermal state of 
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Figure 1: Molecular dipoles orientation (s x (0),s z (0)) at the 
interface boundary depending on the surface interaction pa- 
rameters (ao,/3o). 



bulk water in the model is characterized by the two 
scales: R D = yJC/ (4tt + A), 0.15 nm < R D < 0.25 nm, 
is the size of the strongly correlated molecular cluster, 
and Lt = y/C/A ~ Rd\/£, 1.1 nm < Lt < 1.5 nm 
(Lt ^> Rd), is the size of the largest correlated domain 
within the liquid. 

The interaction of the liquid and an immersed body 
surface is described by fTT] : 



G S = -U/CP* f df(a Q s\ 



•A)*i 



(3) 



where df is the area element of the interface surface F, 
the projections s± = ns and sii = s — sj_n are the normal 
and tangent components of the vector s and n is the unit 
vector normal to the interface surface, directed into wa- 
ter bulk. The dimcnsionlcss phcnomcnological constants, 
ao and /3o> characterize the orientation dependent inter- 
action of the water molecules with the interface. These 
parameters are liquid and surface material specific, and 
should be found from either experimental data or molec- 
ular dynamics simulations. Once all these parameters are 
known, the minimization of total functional G with re- 
spect to the independent variable s (r) leads to the Euler 
equation in the bulk and provides the proper boundary 
condition for the vector field s (r). 

To understand the properties of the liquid interfaces, 
consider the semi-infinite water sample resided in the re- 
gion z > 0, contacting an infinite plane surface z = 0. 
Since there is no external electric field in the system, 
E = Ep = (0,0, — 4irP s z ). The mean field solution 
is obtained via the free energy minimum, using the 
trial function in the form s — (s x ,0,s z )e~ z / LT + 



(s x ,0, s z )e~ z / Rr> , where s x ,' y are the four variable 
parameters. The results of the minimization are repre- 
sented on Figure [T] 

We find three distinctly different types of the solutions, 
depending on the properties of the surface (the parame- 
ters ao and /3rj) ■ If Po is sufficiently large, in region /, 
then the water molecules are polarized along the normal 



to the interface surface, Is 



(1,2), 



0, which corresponds to 
„(i)| 



hydrophilic property. Moreover, \s z | ^ and the polar- 
ization of the liquid extends exponentially into the liquid 
\s z \ ~ er z l R o _ j n region //, the water molecules are po- 
larized along the interface surface, \s z ' '\, \s x | « and 

\s x | ~ e~ z / LT , which is exactly what we expect from a 
hydrophobic surface [13]. In region III, the variational 
solution vanishes and the polarization of the liquid can 
exist only due to the thermal fluctuations. 

Boundaries between the three regions can be found in 
an analytical form, due to the exceptional simplicity of 
the mean field solutions, using the trial function s = 
{s x e- z / LT ,0,s z e-*/ RD ), so that G = SPgy/C ■ (R 2 + 
BR 4 ), 
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where u = y/A, and v 
of G with respect to s a 



= \/47r + A. The minimization 
s XjV shows that the interface is 
hydrophobic (Region I) if /So > v ( a a < u) or Po > 
P(ao) = v + (ao — n) V 'u/v (ao > u)- The wetting en- 
ergy is G/ = -SPl^fCv (p - vf /AB. The hydrophobic 
type of the interface (Region //) corresponds to ao > /i 
and Pa < P(a ), where G n = -SP§VCu (a - u) 2 /AB. 
The fluctuation dominated Region III corresponds to 
ao < u and pa < v, when the mean field Gm = 0. 

Interaction forces between plane surfaces in water for 
Regions / and // can be calculated using the same for- 
malism. Consider first two hydrophilic bodies (e.g. bio- 
logical membranes) with plane surfaces at z = ±d/2, sep- 
arated by the water filled layer of width d. In extreme 
hydrophilic case, ao « and po ^> 1, the mean field 
solution gives s z (±d/2) = ±so, where sq < 1 and the 
minimization of G recovers the experimentally observed 
dependence (fil) with Pi = 2ir PqSq ~ 3-10 10 dyn/cm 2 and 
A = Rd- Similarly in Region II, for d > Rd we obtain 
\s z \ w and P = -P 2 e- d / x , where P 2 ~ P x ~ 2irP 2 and 
A ~ Lt in agreement with the experimental data (5J [3] . 
Therefore, in both cases the interaction force decays ex- 
ponentially with the distance between the planes. Both 
the decay length and the pre-exponential factor depend 
on the properties of the surface material. 

Region III represents a very special case, where the 
mean field polarization vanishes and the energy of the liq- 
uid is determined by thermal fluctuations. The geometry 
dependent part of the free energy leads to the interaction 




mally diverges. To compute the sum, we formally write 



Figure 2: Graphical solution of the characteristic equation for 



between the boundaries exactly in the way electromag- 
netic field fluctuations lead to appearance of the Casimir 
forces [14 . To describe the fluctuations, we use Eqs. 
(|2|-(|3|, keeping only terms ~ s 2 . The free energy of the 



liquid takes a form G 



H 



s ) , where H is a properly 



constructed self-conjugated operator. Diagonalization is 
produced by the decomposition s = ^2 n c n s n over the 
complete set of orthogonal and normalized eigen-mode 
functions {s n }, so that G = J2 n ^n\ c n\ 2 t where k 2 are 
the eigen-numbers of H corresponding to the modes s n 
and enumerated by the index n. Each c n is an indepen- 
dent variable. 

In the case of the two plane surfaces separated by 
the distance d, the liquid is translational invariant along 
the interfaces surfaces and the mode functions can be 
represented in the form s n = (u n (z) ,0,v n (z)) exjp e lpx . 
The solutions are characterized by the set of numbers 
n = (p, P, to), where p is the two-dimensional wave vec- 
tor and P — ±1 is the parity of the function, u n (—z) = 
Pu n (z). Depending on the parity, the mode functions 

(1.2), 



u n (z) and v n (z) are the linear combinations of sing., 

J (1.2) 

and cosq> m 'z, 



where qL = \/k% 



P 



A and g™ = 



47r. In the practically important case d > Rd, all 

the terms containing q, n are small, ~ exp(— d/Ru), and 
as such can not contribute to the interaction force. This 
is because at sufficiently large distances from the surfaces 
only the so-called "force-less", Ep — 0, fluctuations of the 
liquid contribute to thermodynamic functions [9 . The 

wave vector q m = q\n is the solution of the characteristic 
equation: / (X m ) = aX m , where X m — q m dj{2\fC) > 0, 
and a = 2\/C /{ptod). The function / (X) has the period 
7r/2 and is obtained by a periodic shift of tt/2 of the main 
branch of tangent function, tanX for — tt/2 < X < tt/2. 
as shown on Figure [2j 

The equilibrium free energy of the system is 



G eq (d) = -Tin f Ds{v)e- G ' T = ~ f^f] lnk m , 
J O J (2tt) m=1 

and, exactly as in the calculation of Casimir energy, for- 



G eq (d) 
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pdp^2 / 
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Now we can follow [15 and perform the summation us- 
ing the Cauchy's argument principle. After the regu- 
larization, G eq (d) — ► G eq (d) — G eq (oo) , we derive the 
expression: 



G eq (d) 



TS 

AttL 2 t 



tdtln 



(At — tanhyoO tanhyo^ 
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where yo = d/2Lx is the dimensionless distance, and A is 
a material dependent quantity (in Region Iff fi/a Q > 1). 
In the two most important limiting cases, the interaction 
energy takes form 



G eq (d) 



TS 
' 2ixL T d 



exp 



, d> Lr, 



TS 
G e9 (d)«-0.53—^, d<L T , 

corresponding to the attraction. This means that the 
interaction is universal, the dependence on the material 
constants is weak and can only be found at intermediate 
distances d ~ Lt- 

The attraction of hydrophobic bodies in our model has 
an entropic nature, in accordance with earlier predictions 
|16( IT?] . Hydrophobic surfaces order molecules of water 
and the effects of the molecular ordering in pi manifests 
itself in two different ways: from the entropy contribu- 
tion to V^s 2 ), as in our case, or through the Oseen en- 
ergy term, modeling the hydrogen bonding. The latter 
describes the short-range forces and decays at distances 
~ Rd- Hence, for small bodies of sizes < Rjj the hydra- 
tion energy is proportional to the volume |161 \TT\ . The 
longer range interaction between hydrophobic bodies at 
distances ~ Lt originates from the long-range dipole- 
dipole interaction between the molecules and thus re- 
quires a complete model like ((2|-(l3|, which naturally in- 
cludes both distance scales. 

On a side note, we predict that under specific condi- 
tions there could be a special limit. When the liquid 
interfaces fail to polarize water molecules, the fluctua- 
tions of molecular polarization become strong and the 
interaction becomes very weak but attractive. The fluc- 
tuations may also be relevant next to hydrophobic inter- 
faces, where the mean field ordered state of the liquid 
may break in a BKT-like phase transition |111 ITS] , which 
have been observed in molecular dynamics calculations 
for the hydration water layers |19j . 

In summary, we find the polar liquid phenomenology 
(pl-pl earlier proposed |9~HTT] paints a very physically 
rich picture of possible wetting regimes and interactions. 



If used to calculate the interactions between the hy- 
drophilic planes, the model can be considered as the nat- 
ural improvement of ideas [H El [2QH23- Our model gen- 
eralizes the order parameter (the net polarization of the 
water molecules) in a form useful both for hydrophobic 
and hydrophilic interfaces, correctly predicts the sign and 
the distance dependence of the interaction forces depend- 
ing on the properties of the surfaces. We show that the 
interactions originate from the fundamentally non-linear 
and non-local polarizability of the liquid. Our model is 
characterized by two scales, Ru and Ly, instead of a 
single scale ~ R D from [6-8 . The experimental observa- 
tion of both decay lengths in [Tj and [3] together with 



the prediction [TU1 E5] and subsequent experimental ob- 
servation of the ferroelectric features in the bulk liquid 
water near the A— point [121 1251 127] shows both general 
consistency and applicability of the model for realistic 
calculations of macroscopic bodies interactions in water. 
This observation makes our model the minimal continu- 
ous model capable of predicting finer effects depending 
both on the hydrogen-bond network properties and the 
electrostatic interactions of the water molecules. 

The authors are grateful to Prof. V.G. Levadny for 
the fruitful discussions. Quantum Pharmaceuticals sup- 
ported this work. 
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